A MATLAB method of lines template for transport equations

نویسندگان

  • H. S. Lee
  • C. J. Matthews
  • Roger D. Braddock
  • Graham C. Sander
  • F. Gandola
چکیده

Many environmental problems involve diffusion and convection processes, which can be described by partial differential equations (PDEs). This paper will describe the development of a MATLAB template that generates a numerical solution to PDEs using the method of lines. The template will be applied to various unsaturated flow problems within soil physics to demonstrate the versatility of the method. In particular, the template will generate solutions for three cases (1) one-dimensional Richards’ equation for vertical infiltration; (2) coupled one-dimensional Richards’ equation and solute transport equation for horizontal water and contaminant flow; and (3) two-dimensional Richard’s equation for unsaturated flow over a complex geometry. Where possible, the results from the template will be compared against analytical solutions to determine the accuracy of the numerical solution. In addition, the paper will provide a discussion on possible extensions to the template and future directions.  2003 Elsevier Ltd. All rights reserved.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Solving reaction–diffusion equations 10 times faster

The most popular numerical method for solving systems of reaction–diffusion equations continues to be a low order finite–difference scheme coupled with low order Euler time stepping. This paper extends previous 1D work and reports experiments that show that with high–order methods one can speed up such simulations for 2D and 3D problems by factors of 10–100. A short Matlab code (2/3D) that can ...

متن کامل

Sequential Implicit Numerical Scheme for Pollutant and Heat Transport in a Plane-Poiseuille Flow

A sequential implicit numerical scheme is proposed for a system of partial differential equations defining the transport of heat and mass in the channel flow of a variable-viscosity fluid. By adopting the backward difference scheme for time derivative and the central difference scheme for the spatial derivatives, an implicit finite difference scheme is formulated. The variable-coefficient diffu...

متن کامل

تخمین عملکرد خطوط و ریسک عایقی ناشی از ضربات مستقیم صاعقه با استفاده از روش ترکیبی جدید مبتنی بر روش مونت کارلو

This paper presents a new method, based on analytical relations and Monte Carlo simulation for calculating the lightning performance of overhead lines caused by direct strokes. The aim of the work is to reduce the computing time together with retaining enough accuracy to estimate the line lightning performance and performing the risk analysis of power networks. In the proposed method, some modi...

متن کامل

Development of a Computer Model for the Prediction of Littoral Sediment Transport Rate

In this paper, longshore sediment transport in littoral zones is investigated. For investigation of sediment transport in the nearshore zone, the effects of waves, currents and topographical conditions of coast are considered. Linear wave theory was used for the investigation of the wave behaviour. Governing equations of littoral current are continuity and momentum ones. For calculating concent...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • Environmental Modelling and Software

دوره 19  شماره 

صفحات  -

تاریخ انتشار 2004